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Abstract 

All three-loop on-shell QCD Feynman integrals with two masses can be reduced to 
27 master integrals. Here we calculate these master integrals, expanded in e, both 
exactly in the mass ratio and as series in limiting cases. 



1 Introduction 



Massive on-shell Feynman integrals have numerous applications. We consider 
two-leg diagrams, where the external particle has mass M and an on-shell 
momentum p {p^ = M^), in QCD or QEE0. On-shell Feynman integrals with 
a single mass M have been investigated at two p!|2|3] and three [if5|6] loops. 
Starting from two loops, diagrams with loops of massive particles having a 
different mass m appear. Such diagrams are non-trivial functions of the ratio 



(1.1) 



At two loops, there is one generic class of diagranf^ in QCD and QED (Fig. [T]). 
The corresponding Feynman integrals can be reduced to 4 master integrals 

^ These results are also useful in the electroweak theory, or more general field 
theories; but there additional classes of diagrams appear. 

^ We use the following definitions: a class consists of diagrams with identical de- 
nominators, or whose denominators can be made identical by linear substitutions 
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(Fig. [2]), coefficients being rational functions of d = 4 — 2e and x; a reduction 
algorithm has been constructed in Ref. f^f^ . 




— ^--^ 



Fig. 1. The generic class of two-loop on-shell diagrams with two masses. Solid lines 
have mass M, dashed lines mass m, and wavy lines are massless. 
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Fig. 2. The two-loop master integrals. 
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Fig. 3. The generic classes of three- loop on-shell diagrams with two masses. 

The ffist three-loop on-shell calculation with two different masses has been 
done by Laporta and Remiddi [8]. At three loops, there are three generic classes 
of diagrams in QCD (Fig. [3]). Recently, three of us together with M. Stein- 
hauser have demonstrated [9j that the corresponding Feynman integrals can 
be reduced to 27 master integrals (Figs. HHO]) (see also ^lOj). The reduction is 
performed by the C++ program Crusher [llj which implements the Laporta 
algorithm [12] to solve integration by parts identities [13]. This is done exactly 
at general d] simple polynomial operations are done by GiNaC [2], and simpli- 
fication of rational expressions in d and x by Fermat [TB]. On-shell mass and 
wave-function renormalization of a heavy quark (e.g., fe) in QCD with another 
massive flavour (e.g., c) have been calculated at three loops. In the present 
paper, we present details of the calculation of the master integrals used in [9]. 
These integrals can be used in many other three-loop on-shell calculations. 

Our results can also be used for another important kind of diagrams, namely, 
diagrams with two different non-zero masses and any number of external lines 



for their integration momenta; a class is generic if it contains a maximum number 
of lines (denominators); contracting some line(s) in a class, one obtains contracted 
classes. 

^ In this paper, another integral was used as master instead of the last one in Fig. [2l 
It is, of course, easy to reduce this last integral to the master ones of Ref. [7J, solve 
for this master integral, and then express all results via the integrals of Fig. [2j 
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Fig. 6. The master integrals with 4 lines: two-loop sunsets. 
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4.8 4.8a 4.8b 

Fig. 7. The master integrals with 4 lines: three-loop sunsets. 



with small momenta. If we expand in them, these diagrams reduce to vacuum 
integrals with two masses. And such integrals are particular cases of on-shell 
integrals considered in the present paper (the incoming leg and the outgoing 
one are attached to one and the same vertex). Such diagrams reduce to the 
master integrals in Figs. HI [51 Therefore, the master integrals we are investi- 
gating here are also useful for another wide area of applications. 
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We shall use the Minkowski notation. All denominators will contain —iO. Loop 
integrals will be divided by in^^'^ for each loop. For example, 



o 



a 



1 



Via 



(M2 - F - i^Y 



Via] 



-M" 



r(a) 



-m 



(1.2) 
(1.3) 



Powers of denominators are either written near the corresponding lines, or in- 
dicated by dots (if they are some specific small integers). The master integrals 
with 3 lines (Fig. HI) are products of ([L2]), fOll . 
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2 Two-loop sunset diagrams 



The master integrals 4.4-4. 7a (Fig. [6]) are products of two-loop sunset dia- 
grams (Fig. [TUl) and one-loop vacuum bubbles. The first sunset diagram can 
be expressed via F functions for any powers of denominators: 




T{d/2 - ai)T{d/2 - a2)F(ai + 02 - d/2) 
T{ai)T{a2)T{as) 



X 



F(ai + 02 + CI3 — d)r(2{d — ai — 02) 



as, 



T{d — ai — a2)T{3d/2 — ai — 02 — as) 



j^3cZ-2{ai +02+03) 



(2.1) 





3.2 3.3 3.1 

Fig. 10. Two-loop sunset diagrams. 

The non-trivial single-scale master integral (3.3 in Fig. [TOj) was investigated 
in Ref. 121 • It can be written as 



X 



1 - 6e 



(1 -2e)(l -3£)(2-3£) 

7T{l-e)r\l + 2e)T{l + 3e) 
1 F2(l + £)F(l + 4£) 



8_2o-6sr(l + 2e)F(l + 3e) 



F5(l + £) 



+ 2Ae^B^{e) 



M 



2-4e 



(2.2) 



where 



BAe] 



3^2(1 + 2e) 



3F2 



1 i -£ i 

r, C, „ C 



3 3 
2' 2 



7 /F(l -e)F2(l + 2£)F(l + 35) 



24^4 \^ F2(l + £)F(l + 4£) 
7r22-''^F(l + 2£)F(l + 3e) 



3^2 
16Li4 



^V|log^2(log^2 



— TT 



- 1 



13 
180 



(2.3) 
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The expansion of B4^{e) up to O (e^) can be found in p!6]. The hypergeomet- 
ric function can also easily be expanded in e with the help of the package 
HypExp [i7\. 

The two master sunset integrals with two masses (Fig. [2]) have been calculated 
up to in Ref. [18], and up to O (e^) in Ref. [190. Here we obtain 

exact results for them to O {s^)- To this end we use the method of differential 
equations [21]. The class of two-loop sunset diagrams with two masses has two 
master integrals. This leads to a system of coupled differential equations. The 
first one is 



d /I 1 

+ 



dx ^ ^ \x + 1 X — 1 

- 326M' 



X + 1 X — 1 



1 1 



\x X + 1 X — 1 



1 



\x X + 1 X — 1 



/ \ 



\ I 



(2.4) 



where the coefficients have already been expanded in e and only the leading 
term in e for each coefficient is written. As usual the integrals of contracted 
classes are assumed to be known and so are part of the inhomogeneous term of 
the differential equation. If the system is solved order by order in e the integral 
with two dots (3.1b) decouples from (12.41) . The equation for this integral is 

d ^ ^ Ix^ 

-r -< • • — >- = - -< • • >- • 2.5 



Note that here e is set to zero in all coefficients. Of course all integrals which 
appear in ( 12. 4p are also present in this equation, but they only contribute at 
higher orders in e. Looking at the structure of the coefficients in this system 
one can see that it is possible to get solutions for the master integrals in terms 
of harmonic polylogarithms [2^ (this is also true if one writes the coefficients 

^ they were also calculated up to O {e) as series in x up to [20] . 
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to arbitrary order in e) . We have used the package HPL [23] to implement the 
integration of the differential equations. 

After integrating there is one constant per order in e and per integral which 
is not determined. In the case of the integral 3.1 (Fig. [TUl) one constant can 
be determined in the limit x ^ 0, where this integral reduces to 3.2. For the 
integral 3.1b the limit x ^ gives the correct result independently of the 
choice of the constant, as the constant in the integral is proportional to x and 
drops out in this limit. To get this constant we use the limit a; ^ 1, where the 
resulting integral can be reduced to 3.3. 

The result for the master integral 3.1 (Fig. [TU]) up to the finite part in e reads 



1/2 1 



+ ^ (^Ax^ H{0; x) - 3x^ - ^) + 2(a;^ - 1)^ [^(-1, 0; x) - H{1, 0; x) - y 



Ax^{x^ + 2)/7(0, 0; x) + lAx'^ H{0; x)-6x^ -^ + [e] 



(2.6) 



where H denote the harmonic polylogarithms. Up to this order they can of 
course be expressed in terms of logarithms and ordinary polylogarithms, but 
this gets impossible in higher orders in e. We have calculated the diagram 3.1a 
(Fig. [TOj) (which can be reduced to 3.1b) up to order e^, but we refrain from 
presenting the result in a written form in this paper. Instead we refer to the 
web-site |23], where all our results can be found in the form of a Mathematica 
package. Our results agree with [19| . 



We also used another approach. The Mellin-Barnes representation of the mas- 
sive one-loop self-energy is |25p1 



ai 



'^^ iii'^l'^ J [m^ — A;^ — iO]"^ [m^ — (A; + p)2 — iO]' 



r(ai)r(a2) 2Txi 



dzTi—z) 



r(ai + z)T{a2 + z)T{ai + 03 - d/2 + z) -z /o ^\ 

X ^ m •wvwwwv* 2.7 

V{ai + a2 + 2z) ^ ' 



For a simpler derivation, see [26j . 
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where the integration contour is chosen in such a way that poles of V functions 
with +z are to the left of it and of those with —z are to the right. We obtain 



0-2 



J\^d-2a3 ^d-2[ai+a2) 



dzx 



-2z 



X T{-z)T{a^ - d/2 - z) 

r(ai + z)T{a2 + z)T{ai + 02 — d/2 + z)T{d — 03 + 2z) 



X 



T{d — cts + z)T{ai + 02 + 2z) 
This result is a particular case of [7]. 



(2. 



3 Three-loop vacuum bubbles 



The diagrams in the first row of Fig. [5] are vacuum bubbles. The first one, 4.1, 
can be expressed via F functions for any powers of denominators: 




\ _ T{d/2 - ai)T{d/2 - a2)F(ai + 02 - d/2) 
a3, ai^ ^2 ^4 - r(ai)F(a2)F(a3)F(a4)F(d/2) 



X F(ai + a2 + as — (i)F(ai + 02 + 04 — 6?) 



X 



F(ai + a2 + + — 3d/2) 
F(2(ai + 02 - (i) + 03 + tti] 



m 



3(i— 2(01+02+03+04) 



(3.1) 



The non-trivial single-scale master integral (4.2 in Fig. [5]) was investigated in 
Ref. [2]: 



/ ' ^ \ 

/ I 

I I 

1 I 

^ \ / ' 



3(1 -£)(1 -2e)(l -3£)(2-3£) 



F2(l +e)F(l + 4e) 
where B^^e) is given by Eq. (12.31) . 



m 



l-6e 



4-6£ 



(3.2) 
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As in the case of the two-loop sunset diagrams, the vacuum bubble diagrams 
4.3 and 4.3a can be calculated by solving the corresponding differential equa- 
tions in terms of harmonic polylogarithms. We have used the solution of di- 
agram 4.1 (see Eq. fl3.ip ) as initial condition in the limit x ^ 0. Taking the 
limit a; — > 1 we recover the expanded version of Eq. (13.21) . For example, 




+ 4x^ + 1 

3^3 



- ^ (2x^{%^ + 2)if(0; x) - ^(7a;^ + 32^' + 7)) 
+ ^ [^x^{^x^ + 2)/7(0, 0; x) - x^{7x^ + 16)H{0; x) 



+ — (5x^ + 32x2 + 5)) +C(1) 



(3.3) 



The analytic results expanded up to O (e^) of all vacuum- diagrams depicted 
in Fig. ([5]) can be obtained on the aforementioned webpage. 

For the vacuum bubble with two masses we obtain, using the Mellin-Barnes 
representation (12.71) . 




■j^2{d-a3-a4)^d-2{ai+a2) 



dzx 



^2z 



X 



r(ai)r(a2)r(a3)r(a4)r(rf/2) 2m 

r(— 2;)r(a3 — d/2 — z)T{ai — d/2 — z)T{az + — d — z 

r(a3 + — d — 2z) 
T{d/2 + z)T{ai + z)T{a2 + z)T{ai + aa - d/2 + z) 
Tiai + a2 + 2z) 



(3.4) 



(this is a particular case of a more general result recently derived in [27]). This 
diagram is, of course, symmetric with respect to ai^2 ^ '^3,45 M ^ m, though 
this is not evident in (13.41) . 



4 Three-loop sunset diagrams 



Here we consider three-loop on-shell sunsets with masses M, m, m, (Fig. [7]). 
The first of them has been calculated up to O (e^) in Ref. [28], though only 
the O (e^) term is presented in the paper. 
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Doing a naive reduction to master integrals of all three-loop on-shell sunsets 
with masses M, m, m, one finds four "master integrals" . As already noticed 
in [28j, one of these integrals decouples completely when setting up a higher 
order differential equation for the integral 4.8 (Fig. [7j). We have performed 
a reduction for all classes of the different master integrals presented in this 
paper. It turns out that in the course of the reduction of the class to which 
the integral 5.1 (Fig. [8]) belongs, an equation is generated which contains no 
five-line master integrals and connects the integral 4.1 (Fig. [5]) and one of the 
integrals from the class discussed in this section. So it gets manifest that one 
of the four integrals found by the naive reduction is reducible with integration 
by part identities. 



For example, for 4.8 we obtain 



1 



3^3 



+ — -2x\x^ + l)H{0; x) + -{7x^ + 15x^ - 1) 

\ 6 

+ - (^4x^(3x2 + l)H{0, 0; x) - 7x^{x^ + l)i/(0; x) 

+ T^(75a;^ + 213x2 - 35)) + O (1) 
3o / 



(4.1) 



The analytic results for all three master integrals in Fig. [7| up to order e are 
presented on the aforementioned webpage. 



Using the Mellin-Barnes representation (12. 7p . we obtain 



^ a2 ^ 

/ _ - - - N 



V.,,,^^^^,^^ r(ai)r(a2)r(a3)r(a4) 
as 

X — f dz x'^'r{a3 - d/2 - z)r{a3 + - d - z)r{d/2 + z) 
2m J 

— joo 

^ r(ai + z)T{a2 + z)T{ai + a2 - d/2 + z)T{2{d - as) - a4 + 2z) 
T{d — as + z)T{3d/2 — — + z)r(ai + a2 + 2z) 
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5 Master integrals with 5 lines 



We have calculated the master integrals 5.1 and 5.1a (Fig. [8]) analytically up 
to O (e) with the differential equations method. The integrals 5.4 and 5.4a are 
obtained with the same method up to O (e^). The integrals 5.2, 5.2a, 5.3, 5.3a 
are obtained up to O {e~^) with that method. 

For the two classes 5.1 and 5.4 the differential equations have a structure 
with which it is possible to integrate them with the help of harmonic poly- 
logarithms. So in principle, if the initial conditions are known for the corre- 
sponding integrals, solutions can be obtained to arbitrary order in e. This is 
different for the classes 5.2 and 5.3. Here the pole structures of the differential 
equations do not only contain poles of the form 1/x and 1/(1 ± x), but in 
addition the poles 1/(1 ± 2x) (for 5.2) and 1/(1 ± x/2) (for 5.3). It is not 
possible to integrate the differential equations in terms of the usual harmonic 
polylogarithms with three weight functions. Nevertheless we were able to in- 
tegrate all equations for these integrals up to order e~^. The 5° parts of the 
master integrals 5.3 and 5.3a were obtained analytically using Mellin-Barnes 
representation, see below. 

To get a result in higher orders in e for the integrals of the classes 5.2 and 5.3 
we calculated the integrals in an expansion around x = 0. With the help of 
the ansatz 

^ Cijfc e V log'^ a; (5.1) 

for the master integrals, the differential equations can be expanded in e and 
x. As a result the differential equations reduce to algebraic equations for the 
coefficients Cijk- In every order in e there is one constant Cijk which can not be 
determined with this procedure. 

A problem which can arise when determining these remaining constants is the 
fact that it may happen that the initial condition for some integrals is fulfilled 
for the corresponding constant being arbitrary. In the case of the class 5.3 this 
problem occurs when one chooses the two integrals depicted in the third line 
of Fig. [H] as master integrals. Sometimes one can determine the constants via 
analytic considerations. In our case this was not possible. The problem can be 
solved by switching to another master integral basis and replacing the integral 
with the dot with the integral 5.3b: 
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in 



M2 - (/3 - i,y 



where = denotes the on-shell external momentum. For this particu- 
lar choice of master integrals it is possible to determine all constants in the 
ansatz (15. ip using the corresponding initial conditions at x = 0. 

For the class 5.2 no basis of master integrals exists, where all constants can be 
determined via the initial conditions. In every order in e the coefficient of the 
x^-term in the x-expansion of the integral 5.2 remains undetermined. To get 
this constant we calculated the integral by means of the method of regions [21] 
up to the third order in x. After having determined the constant in this way 
all other orders in x are fixed. 



On the website mentioned earlier we present the analytic results for the in- 
tegrals depicted in Fig. ([8]). We also give expansions up to O (x^^) for the 
integrals depicted in the second and third line of the figure up to O (e). 

For the diagrams 5.1 and 5.1a in Fig. [HI we can obtain a one- fold Mellin-Barnes 
representation, using (12.71) : 



14 a^. ^ \ 



as 



+ioo 



T{d/2 - a3)T{d/2 - a4)r(a3 + - d/2) 1 1 , -2z 



r(ai)r(a2)r(a3)r((24)r(a5)r((i — — a^) 27ri 

X r(— z)r(a3 + + — d — z)T{ai + z)T{a2 + z) 
r(ai + a2 — d/2 + z)V{2{d — as — a4) — as + 2z) 



X 



r(3(i/2 — as — a4 — as + z)T{ai + a2 + 2z) 



(5.3) 



For the diagrams 5.2, 5.2a, 5.3 and 5.3a in Fig. [SI we have to use (12. 7p twice. 
Therefore, we obtain two-fold Mellin-Barnes representations: 
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r(ai)r(a2)r(a3)r(a4)r(a5) 



Hod +100 



X 



(27rz) 



dzi / rf;z2a;-^("^+"^) 



-loo —too 



X 



X 



X 



r(-zi)r(ai + 2i)r(a2 + -2i)r(ai + 02 - (i/2 + 2;i) 

r(ai + 02 + 22:1) 
r(-Z2)r(a3 + Z2)r(a4 + Z2)T{a3 + - d/2 + Z2] 

r(a3 + a4 + 22:2) 
r(a5 - d/2 - zi - Z2)T{d - as + 2{zi + Z2)) 
T{d - as + ^1 + Z2) 



j^^2{ii-a3 -04-05)^^-2(01 +02) 




r(ai)r(a2)r(a3)r(a4)r(as) 



X 



X 



X 



X 



(27ri 



-lOO — JOO 



r(-zi)r(ai + Zi)r(a2 + -2i)r(ai + a2 - d/2 + zi] 

r(ai + a2 + 22;i) 
r(-Z2)r(a3 + Z2)r(a4 + ^2)r(a3 + a4 - d/2 + Z2) 

r(a3 + a4 + 22:2) 
r(a5 - d/2 - zi - Z2)T{d - as + 2{zi + ^2)) 
T{d - as + 2:1 + Z2) 



(5.4) 



(5.5) 



We were able to get the e -part of the integrals 5.3 and 5.3a analytically 
from (15. 5p . Since we have a reduction of this family to two master integrals 
we could evaluate any pair of linearly independent integrals. We have cho- 
sen Feynman integrals with the indices (1, 2, 1, 2, 1) and (1, 1, 1, 1, 2) (see the 
enumeration in Eq. (15.51) ). We evaluated it using (15. 5p and resolving singu- 
larities in e in the corresponding two MB representations [Sm31f26j with the 
help of the Mathematica packages MB.m [32] and MBresolve.m [33]. In both 
cases, to evaluate the £°-part of these integrals, we needed to evaluate at most 
two-fold finite MB integrals. One of the integrations was done by corollary of 
Barnes lemmas (see Appendix D of [20]) implemented in MB .m. At the last step, 
one-fold MB integrals were evaluated by closing the integration contour and 
summing up series. The expansions in x agree with those obtained from differ- 
ential equations. The results were also checked numerically using the program 
FIESTA which implements sector decomposition, where various algorithms 
of sector decomposition developed in Refs. [S51l3tif34j are implemented. 
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For the calculations presented in [9] we also needed the coefficients of the 
five line master integrals. As discussed above we were able to calculate all of 
them as series in x but not all could be evaluated analytically. Let us now 
describe how we can get a numerical solution for these integrals, including the 
coefficient. We also used this method as a numerical check of our analytical 
results. 

We choose integral 5.3 (Fig. [8]) with all indices equal to 1 as an example. The 
starting point is the MB representation, Eq. (15.51) . of this integral. Inserting 
d = 4 — 26 and a,- = 1 we obtain 




M' 



2-6e 



+ioo 

2e~2zi 



j dzi j dz2X 



-lOO ~ioo 



X 



(27ri)2 

T{-z^)^\l + z^)^{e + z^)r{-Z2)T\l + Z2)T{e + z^) 



r(2 + 2^i)r(2 + 2^2)r(3 -2e + zi + Z2) 
X r(-l +e-zi- Z2)r(3 -2e + 2zi + 2^2) • (5.6) 

The next step is the resolution of the e singularities using the algorithm de- 
scribed in [30|33] . Our integral is decomposed into five analytical expressions, 
two one-dimensional MB integrals and one two-dimensional one, which con- 
tains the same integrand as Eq. (15.61) but with modified integration contours. 
The integrand expanded in e reads 



(x2)-^i 



r (zi + 1) r(^i + ^2 + 3)r(222 + 2) 

X 2^^^+^ r(-^i) V{z^) T{zi + 1) r(-^i -Z2-I) T{-Z2) 

X T{Z2) T\Z2 + 1) r (^Zi + -22 + ^) T{ZI +Z2 + 2) 

X {-1 +£[21oga; + 21og2 - ^{zi) -ip{-zi - Z2 - I) - ip{z2) 
+ ij[zi + Z2 + ^ + i){zi + Z2 + 2)- 2ij{zi + ^2 + 3)] }. 



(5.7) 



To evaluate the integral we close the integration contours and sum up the 
residues. Since x"^ appears with the power —zi and we are interested in an 
expansion for small x we have to close the contour for zi to the left. We are 
free in the way how we close the contour for Z2 and choose the right side, 
considering the F functions F(— 2:2) and F(l — zi — Z2). The integral is then 
given as the sum of the residues in the points Z2 = n and Z2 = —zi -t- 1 -|- n for 
positive integers n. The expressions for these residues read (we write only the 
part here): 
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22n+l ( 2\-zi 

A,(n) = -Res(F, = n) = , ^ 

riTT [n + Zi + 1) 

T\n + 1) r(i - ^i) r(-n - ^1 + 1) r(-^i) r^(^i) r(^i + 1) 

T{2n + 2) r (^1 + I) T{n + Zi + 1) T{n + + 3) 
xr(n + zi) r(n + zi + ^^ r{n + ^1 + 2) sm(7r(n + z^)) , (5.8) 

A2{n) = -Res(F, Z2 = -Zi + l+ n) - ^ ' 



22^1+1 n(n -Zi-1) 

r (n + i) r(n + 1) r(i - T\n - zi) r(-^i) r2(2i) t{zi + 1) 
r(n)r(n + 2)r(2n - 2zi)r (^i + 1) 



(5.9) 



Now we take the residues in the points Zi = —n — | — m and zi = —m with 
positive integers m. We have to choose m > or m > 1, depending on how 
the integration contour passes the poles of the V functions. 

Thus we get a two-fold sum, 

^^/,(m,n), (5.10) 



the summands fi consisting of rational functions of m, n and F, if) and trigono- 
metric functions whose arguments are linear combinations of m, n. 

To check our analytical series expansion we proceed as follows. Some of the 
summands are proportional to (a;^)™'"'""' and some to (x^)™. In the first case 
we only need to sum up a finite number of terms to get an expansion to a 
given order in x. In the second case the sum is an expansion in x where the 
coefficients are infinite sums in n. We evaluate these coefficients to a finite 
order given by the required numerical precision. 

The one-fold sums contributing to the integral 5.3 are special cases of the 
described procedure. However, there is one partial sum which does not depend 
on X and shows quite a bad convergence behaviour. In this case we use the 
method of nonlinear sequence transformations (see Refs. [37|l38] and references 
therein) to improve the convergence. 

As a result we finally get an expansion for the integral 5.3 in powers of x with 
numerical coefficients. Up to the order x^ it reads 



15 



\ 

\ 



^2-6£g-37B£|i_(_o.66667 - 0.33333 x^) 



+ 



+ 



+ 



1 



-3.33333 + (-2.00000 + 2.00000 logx) x^' 

13.60147 + (-0.57606 + 12.00000 logx - 2.00000 log^ a;) 
+ (-5.03987 + 3.00000 logx - 2.00000 log^ x) x'^ 
-52.02282 

+ (31.35527 + 28.61586 logx - 12.00000 log^ x + 1.33333 log^ x) x^ 

- 52.63789 

+ (13.71171 - 0.92026 \ogx - 5.00000 log^ x + 4.00000 log^ x) x^ 
+ 10.52758 x^" 

-153.61196 + (149.09532 + 91.43326 logx 

- 28.61586 log^ X + 8.00000 log^ x - 0.66667 log^ x) x^ 
+ (-94.12475 + 210.55156 logx) x^ 

+ (-106.98907 + 20.56599 log x - 6.01454 log^ x + 6.00000 log^ x 



- 4.66667 log^ x)x^ + (-18.37249 - 42.11031 logx) 



(5.11) 



where we arbitrarily display five digits. We determined the coefficients with a 
relative error of at least 10~^. The result numerically agrees with the analytical 
expansion obtained with the differential equation method. 



To find numerical solutions for the master integrals we insert fixed values for 
X in the sum (15.101) and sum it up numerically. As a further check we use the 
numerical integration routine of MB. The summation procedure takes more 
computing time but its advantage is that the errors are smaller than those we 
get from the numerical integration. 



We used this method to find numerical expansions and also numerical results 
for fixed values of x for the master integrals 5.2, 5.2a, 5.3 and 5.3a. Let us 
finally mention the differences in the calculation of these. Diagram 5.2 has two 
loops with the light quark mass. Therefore the MB integrand is proportional 
to x^^^^'^ and we have no freedom to choose how we close the integration 
contours. On the other hand, all summands contain the parameter x and thus 
the calculation is faster. Integrals 5.2a and 5.3a contain an additional power 
on one denominator. This changes the arguments of the F functions and we 
have less MB integrals to consider. Their complexity is comparable to that of 
the undotted integrals. 
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6 Master integrals with 6 lines 



Contrary to the integrals with five lines, the ones depicted in Fig. [9] pose 
no difficulties when solving them with the help of the differential equation 
method. The integrals where some parts of the e-expansion are only known as 
expansions in x (see the previous Section) do not appear as contracted classes 
in the differential equations. Furthermore, the poles of the coefficients in the 
differential equations are only of the form 1/x and 1/(1 ±x). We can therefore 
get closed solutions of the equations in terms of harmonic polylogarithms. 

To obtain the integration constants we use the boundary condition at x = 1. 
The corresponding values can be found in [5] up to C [e). We obtain for 6.1, 
6.2, 6.3 




= M-^^r=^(l + £) |^^ + C(1)J . (6.1) 

— 

The results for the master integrals with six lines up to O {e) can be obtained 
on the website. 

We were able to derive a one-fold MB representation for the master integral 
6.1 in Fig. [91 Using (12. 7p we get a two- loop on-shell integral with a single non- 
integer index —z. It can be reduced to trivial ones by integration by parts, 
and we obtain 
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]^^-4£^4-2(ai+a2+e) ]_ 

{l-2e)T{ai)T{a2)2^i J 



+ioo 



dzx' 



-22 



T{-z)T{ai + z)T{a2 + z)T{ai + a2 - 2 + e + z) 



Tie 



{l-2e + z)T{ai + a2 + 2z) 
T{e - z)T{2 -2e + 2z) V{e)V'^{\ - e) V{2e - z)V(2 - 4e + 2z) 



+ r(i 



r(2-2e + 2) r(l-2£) r(2-3e + 2) 

r(e - z)Y{2e - z)Y{\ - e + z)Y{2 - 4^ + 2z)' 
^'^ r(-z)r(2 - 2£ + ;2)r(2 - 3e + 2) _ 



(6.2) 



7 Conclusion 



The status of our knowledge of the master integrals is summarized in the 
Tables [THll Here DE means that a result is known analytically, exactly in 
X, and has been obtained by the method of differential equations; MB — the 
same, but using Mellin-Barnes representation; x — expansion of the result (up 
to x^^) is known analytically, from differential equations (and a numerical value 
can be calculated for any x by integrating the MB representation numerically); 
a reference means that a result has been obtained by others. 

Table 1 

Two- loop sunset master integrals (Sect.[2|) 







3.3 


3.1, 3.1a 




DE 


DE 


£3 


DE 


DE 


£4 


m 


IE] 


£5 


m 


m 


£6 


m 





These master integrals can be used for calculating any three-loop diagrams 
with two legs of a massive particle (with mass M) both of which are on the 
mass shell, and there is a loop of another massive particle (with mass m). 
They can be downloaded from [24J in the form of Mathematica files. The first 
physical application, namely the influence of 7^ on the on-shell mass 
of the h quark and its wave- function renormalization constant, has already 
been published p]. We plan to use these results in some further problems. 
A subset of master integrals considered here is also necessary for another 
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Table 2 

Three-loop vacuum and sunset master integrals (Sects. [3l |3|) 





4.2 


4.3, 4.3a 4.8 


f 

, 4.8a, 4.8b 




DE 


DE 


DE 


£2 


DE 


DE 


DE 


£3 


m 




[28] 


£4 


m 






£5 


m 







Table 3 



Three-loop master integrals with 5 lines (Sect. [SJ 

















— — 




— 7* — 




5.1, 5.1a 


5.2, 5.2a 


5.3, 5.3a 


5.4, 5.4a 


£-3 


DE 


DE 


DE 


DE 


£-2 


DE 


DE 


DE 


DE 


£-1 


DE 


DE 


DE 


DE 


1 


DE 


X 


MB 


DE 


£ 


DE 


X 


X 


DE 


£2 








DE 



Table 4 



Three-loop master integrals with 6 lines (Sect. [6]) 





A-:. Ak Ak 




6.1 


6.2 


6.3 


£-3 


DE 


DE 


DE 


£ 


DE 


DE 


DE 



kind of problems — diagrams with any number of external lines having small 
momenta and containing loops of two different massive particles. There are 
several interesting physical problems which involve such diagrams. 

We are grateful to M. Steinhauser for collaboration and to P. Mastrolia 
and E. Remiddi for providing the complete results of Refs. [I9l[28] . 
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